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ABSTRACT 

We explore the evolution of the morphology density relation using the COSMOS-ACS survey and 
previous cluster studies. The Gini parameter measured in a Petrosian aperture is found to be an 
effective way of selecting early- type galaxies free from systematic effects with redshift. We find galaxies 
are transformed from late (spiral and irregular) to early (E+SO) type galaxies more rapidly in dense 
than sparse regions. At a given density, the early- type fraction grows constantly with cosmic time, but 
the growth rate increases with density as a power law of index 0.29 ±0.02. However, at densities below 
100 galaxies per Mpc 2 no evolution is found at z > 0.4. In contrast the star-formation-density relation 
shows strong evolution at all densities and redshifts, suggesting different physical mechanisms are 
responsible for the morphological and star formation transformation. We show photometric redshifts 
can measure local galaxy environment, but the present results are limited by photometric redshift 
error to densities above £ = 3 galaxies per Mpc 2 . 

Subject headings: galaxies: evolution — galaxies: formation — cosmology: observations — cosmology: 
large scale structure — galaxies: structure — galaxies: clusters 



1. INTRODUCTION 

The correlation between galaxy properties such as color 
and morphology with galactic environment, specifically 
the lo cal density of galaxies, was first noted by iHubbld 
(1926). Later, low redshift surveys showed a high per- 
centage of early-type (E+S0) and redder galaxies in clus- 
ters compared with the fiel d which was quantified as the 
morphology-density (T-S1 dDres slcr 1980; Dressie r" et al. | 
19971 ITreu et al.ll2003l : iSmith et alJl2005t Pos tman et all 
20051) and star formation rate-density (SFR-S): (lOemleri 
19MlMehiTck fc Sargendll977l : llutcher k Oemlerlll984l : 
Kauffmann et al.l2004D relations. Spectroscopic redshifts 
from the Sloan Digital Sky Survey (SDSS) enabled an 
extension of these co rrelations over five orders of magni- 
tude in £ at z < 0.1 (|Goto et al.ll2003l : iKauffmann et all 
2004). Such studies clearly indicate the present proper- 
ties of galaxies have been strongly influenced by galaxy 
environment, presumably by galaxy interactions. 
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try a nd spectroscopy (ICowie et al 1 
1994 iDickinson eFatl 120031 



Most investigations of galaxy evolution have fo- 
cused on the star formation rates and stellar 
masses as determined from inte grated photome- 

1991 iLillv et al l 
Glazebrook et alj 1 2004 : 
Bauer et al.ll2005t iDrorv et al.ll2005l : iFeulner et~atll2005 ; 
Juneau et al.| [2005) . In these studies star formation ac- 
tivity is seen to progress to less massive galaxies as the 
unive rse ages - a phen omenon often called "downsiz- 
ing" (jCowie et al.lll999f l. Nevertheless, the process of 
star fo rmation in galaxi es is complex and poorly under- 
stood (|Somervillel 12005) . On the other hand, the mor- 
phology of a galaxy is linked to the angular momentum 
distribution which is simpler to quantify, and hence a 
valuable alternative to star formation. The presence of 
a disk clearly indicates a dynamically-cold stellar pop- 
ulation that has not been significantly disturbed, while 
it is likely that spheroidal systems have been heated by 
interactions, merging or other mechanisms. 

At redshifts above z > 0.3 ground-based imaging 
generally has inadequate resolution for morphological 
surveys since the ty pical size of a galaxy is <0.75" 
(|Ferguson et al.ll2004|) . so HST data is required. The 
small field of view available on HST limited previous 
studies of the T-E relation at z > 0.3 to a few previously 
known galaxy clusters and field samples (|Dressler et al.l 
119971: ITreu et al.ll2003t iPostman et al.ll2005fh Combining 
these studies. ISmith et al.l (|2005l ) found high density and 
low density regions evolved differently. They postulated 
different formation times were possibly responsible for 
the observed effect, but were limited by the look-back 
time and density resolution of the available data. The 
COSMOS ACS data is the first HST imaging survey able 
to study statistically significant samples of galaxies over 
a range of environments and redshifts with high resolu- 
tion morphologies. 

In this paper, we use the COSMOS data, combined 
with photometric redshifts, to quantify the evolution 
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of the T-E relation as a function of redshift over the 
range z = 0.3 to 1.2. These measurement s are combined 
with previous studi e s in the literature (Dressier! 119801: 
Dressier et ail Il997t iPostman et~aH 120051: ISmith et all 



TABLE 1 

Number of Objects in Redshift Bins 



20050 to provide a larger dynamic range in density and 
span in redshift than possible with COSMOS data alone. 
In Section 2 we describe the imaging data and photomet- 
ric redshifts (|Mobasher et al.l [20071 ) used for this study 
and the selection of our sample. The morphological pa- 
rameters used to characterize the galaxies are discussed 
and tested in Section 3. Details of estimating density 
with photometric redshifts and correspondingly poor line 
of sight distances, are discussed in Section 4. Finally, de- 
pendence of morphology on density and redshift is pre- 
sented and discussed in Sections 6 and 7 respectively. 
Throughout this paper we use a standard cosmology with 
fl v ~ 0.7, Q m = 0.3, and H Q = 75. All literature values 
are converted to this cosmology unless otherwise noted. 

2. IMAGING AND REDSHIFTS 

The Hubble Space Telescope Advanced Camera for 
Surveys (HST-A CS) images, taken as a part of the 
COSMOS survey ([Scoville et aU2007bl:lKoekemore et all 
12001 are used to measure morphologies. These im- 
ages cover an area of ~1.8 square degrees in ^81414^ 
and 81 square arc-minutes in FA75W with single orbit 
exposures. The completeness is 50% in F814W for a 
galax y 0.5" in diameter w ith F814W ab — 26 magni- 
tude (|Scoville et al.ll2007bD . The point source depth is 
~ 2 mag deeper but this is not relevant to morphological 
studies. The median image quality is 0.05" and 0.08" 
(FWHM) in F475W and F814W respectively . 

Photometric redshifts ar e determined from multi-band 
ground-based photome try (|Capakil2006D as described in 
iMobasher et al.l (|2007l ). Our present investigation is fo- 
cused on redshifts z = 0.3 to 1.2 for which the photo- 
metric redshifts have an accuracy <J Z /(1 + z) = 0.03 for 
galaxies with Iab < 24 mag. The redshift accuracy was 
determined directly by comparison of the photometric 
redshifts with spectroscopic redshifts for over 900 galax- 
ies at z < 1.2 with l AB < 24 in the COSMOS field. 

A magnitude cutoff of F814W^ < 24 was adopted based 
on experimentation in making morphological measure- 
ments. Below this magnitude, the low surface bright- 
ness also leads to incompleteness for a typically sized 
z = 1 galaxy which i s 0.75"in size ([Ferguson et al.l[2004l : 
IScoville et ai1l2007bft . 

Stars are removed from the objec t catalog using the 
SExtractor iBertin fc Arrioutsl JT996) CLASS.STAR pa- 
rameter measured on the ACS F814W image. Objects 
with CLASS.STAR > 0.9 are considered to be stars. 

A total of 120,187 galaxies meet our magnitude cut 
of ^814^ < 24 and fall outside of the masked regions 
on the images. Of these, 32,958 are meet the criteria to 
be used in our morphology-density analysis (see Section 4 
and 6). TableQ]gives the number of objects as a function 
of redshift . 

3. MORPHOLOGICAL CLASSIFICATION 

Galaxies exhibit a range of morphologies which is dif- 
ficult to quantify automatically. So, classification by eye 
is often used to test the efficacy of automated classifiers. 
At high redshifts two additional problems arise : surface 
brightness dimming (oc (1 + z) 4 ) reduces the visibility 



Redshift 
Range 



F8UW < 24 



M v < -21.2+ 
Look Back Time 



0.2 < z < 0.4 
0.4 < z < 0.6 
0.6 < z < 0.8 
0.8 < 2 < 1.0 
1.0 < z < 1.2 
1.2 < z < 1.4 



18250 
17419 
27018 
19062 
13075 
7975 



2154 
2817 
6193 
7829 
7799 
6166 







1 






m # 



Fig. 1. — Six galaxies with early- type visual classifications and 
late type automated classification are shown. The grey scale images 
are the .F814W images used for the visual classification with a 
square root stretch. The color images have the same stretch and 
cuts, but are made from the HST-ACS F475W and F814W images 
which were adaptivel y smoothed using the method described in 
Scovillc et al. (2007a). Notice the low surface brightness features 
obvious in the color images, but hidden in the grey scale images. 



of disk identifying features such as arms and bars and 
redshifting (or band shifting) means that the observed 
visible bands increasingly sample the rest frame ultravi- 
olet (UV). The UV is unlikely to provide a reliable clas- 
sification of the major galactic stellar components since 
it is highly biased in favor of the youngest star forming 
regions and areas of low extinction. These redshift de- 
pendent effects introduce systcmatics which are difficult 
to separate from the desired evolutionary trends. 

Most o f the automated mor p hologi c al classification 
schem es ([Abraham et al.l 11994 Il996t IConselice et al.l 
120031 ) have focussed on separating elliptical, disk and ir- 
regular systems; these galaxies exhibit a wide variety of 
structure and their appearances are very dependent on 
viewing angle and other projection effects. On the other 
hand, dynamically hot systems (E+S0 galaxies) are gen- 
erally more centrally concentrated (therefore less suscep- 
tible to surface brightness dimming) and more spheri- 
cal (therefore less subject to projection effects). Thus, 
many difficulties in the morphological classification dis- 
cussed above are minimized by selecting just the early- 
type (E+S0) population. Specifically, independent of the 
method, we classify galaxies with respect to whether they 
are or are not early-type and do not differentiate within 
the late type population. This binary classification is en- 
tirely adequate for studying the T-E relation so long as 
any systematic biases do no vary with redshift or density. 
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Fig. 2. — The difference between logio (asymetry) measured in 
F475W and F81AW is shown with redshift. Notice the systematic 
trends with redshift as the two bands go from sampling the rest 
frame optical to the rest frame UV. Clumps of active star formation 
are more visible in the rest frame UV, increasing the measured 
asymmetry. This transition happens at z ~ 0.2 for FA75W and 
z ~ 1.0 for F814W. These transitions are indicated with dotted 
lines. 



Fig. 3. — The difference between logio(Gini) measured in FA75W 
and .F814W is shown with redshift. Notice there are no systematic 
trends with redshift as the two bands go from sampling the rest 
frame optical to the rest frame UV. The transition from optical to 
UV light happens at z ~ 0.2 for FA75W and z ~ 1.0 for F814W. 
These transitions arc indicated with dotted lines. 



Surface brightness has proven particularly problematic 
for previous morphological studies (R. Ellis & R. Abra- 
ham private communication). Both visual and auto- 
mated classification schemes miss low surface brightness 
disks, resulting in systematic effects with magnitude and 
redshift. These systematic effects can be mitigated by 
defining the edge of an object in a way which is inde- 
pendent of signal-to- noise per pixel. Figure Q] illustrates 
one such method for eyeball classifications. The origi- 
nal HST-ACS F814W images are shown together with 
adaptively smoothed images for several galaxies with 
faint d isks, all of which wer e visually classified as early- 
types ()Scoville et al.|[2007al ). The stretch and scale fac- 
tors are identical for all images. The adaptive smooth- 
ing gives equal signal-to-noise per resolution element but 
with varying spatial resolution dependent on the local 
signal to noise ratio - i.e. higher resolution in brighter 
areas. This enhances the visibility of low surface bright- 
ness, extended features. In the processed images, spiral 
structure becomes visible in all four galaxies although it 
was difficult to see visually in the original images. 

Petrosian apertures provide object limits (edges) inde- 
pendent of signal-t o-noise per pixel for automated classi- 
fication methods (|Lotz et al.l |2004) . We use a 'Quasi- 
Petrosian' aperture to minimize the effects of surface 
brightness dimming. This was constructed using a new 
algorithm that is intended to work for galaxies of arbi- 
trary shape, and which has more graceful convergence 
properties than the usual formulation of Petrosian aper- 
tures. (Ordinary Petrosian indices are based on circular 
apertures, and are not guaranteed to converge) . The 
details of our procedure are given in ([Abraham et al.1 
l2007b| ) and only an outline is given here. The first step 
is to use SExtractor to isolate the galaxy from the sky. 
The flux values of the pixels in the galaxy are then sorted 
in decreasing order to construct of curve of sorted flux 
values. This curve is then summed over to construct a 
curve of cumulative flux values, which is then multiplied 



Comparison of Qualitative and Quantitative Morphologies 
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Fig. 4. — Morphological parameters measured from the F814W 
HST images are plotted for objects with visual morphologi- 
cal classification. The Gini parameter and rotational asymme- 
try measured in a Petrosian aperture are plotted for objects 
with visual classifications. Dividing lines are drawn between 
regions of predominantly irregular, spiral, and elliptical types. 
The dashed line between irregular and spiral galaxies is de- 
fined as login (Asymmetry)= 2.353*login (Gini)+0.353, while the 
solid line between spiral and early-type galaxies is defined as 
login (Asymmetry)= 5.500*logm (Gini)+0.825. The dotted line is 
a cut at login (Gini)=-0. 35, used to select early-type galaxies in 
this paper. This selection gives an equivalent fraction of early- type 
galaxies as the Gini-Asymetry selection at 0.2 < z < 0.4. Stars 
and very compact objects appear in the upper right of this plot, 
but are removed in our later analysis. 

by a scale factor r\. The flux value where the scaled 
cumulative distribution intersects the sorted pixel value 
distribution defines a critical flux value. Only those pix- 
els with fluxes greater than this critical value are used in 
calculating the 'Quasi-Petrosian' coefficients. 

Rotational asymmetry and Gini are two commonly 
used parameters for automated morphological classifica - 
tion ([Abraham et alj 11994 [l99l IConselice et al.ll2003j) . 
The asymmetry coefficient is calculated from the differ- 
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Fig. 5. — A magnified region of Figure [4] around the division of early and late-type galaxies is shown. The points are replaces with color 
images of the classified objects. The white line corresponds to the division of early and late type galaxies in Gini and asymmetry, while 
the yellow line is our cut at logio (Gini)=-0.35. Notice the clear division of early and late type galaxies by our cut in Gini (yellow line). 



ence between the galaxy image and it's image rotated 
180° about it's central peak. The Gini parameter mea- 
sures segregation of light into bright and faint pixels and 
is strongly correlated with the concentration of galaxy 
light, but the Gini parameter is a more robust indicator 
of gala xy morphology than th e usual concentration coef- 
ficient (jAbraham et al.l l2003f ) . In particular Gini is less 
sensitive to surface brightness effects and does not require 
a well defined centroid. Specific details of our Gini and 
asymm etry measurements are given in lAbraham et al.1 
(20073) • 

The F814W HST data which covers the whole field 
will sample a different rest frame wavelength range at 
each redshift. So any morphological classification scheme 
used to study galaxy evolution must be independent of 
the rest frame wavelength. We can test for systematic 
variations due to band shifting in our data by using the 



the central 81 square arc minutes of COSMOS with dual 
band coverage (F814W and F475W). Figure [2] shows the 
difference between asymmetry measured in F475W and 
F814W bands with redshift. A shift between the classi- 
fications derived in the two bands is observed at z ~ 0.2 
where the median wavelength of the F475W band moves 
into the rest frame ultraviolet (UV) light while F814W 
still samples the rest frame visible. A shift back to con- 
sistent classification occurs at z ~ 1.0 where both bands 
sample the rest frame UV light. As noted earlier, the rest 
frame UV selects regions of active star formation which 
tend to be irregular/clumpy; as a result, the asymme- 
try is significantly higher in the rest frame UV than in 
the rest frame optical. Figure [3] shows the difference be- 
tween Gini measured in F475W and F814W bands with 
redshift. The scatter is much smaller, and no systematic 
trends are observed with redshift. 
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Fig. 6. — Comparison of real-space and projected density estimates for a simulated catalog with structures of known density. The left 
panel (a) includes no redshift error, the right panel (b) includes a redshift error of <r z /(l + z) = 0.03. The dashed line in panel b indicates 
the minimum measurable density, £ = 3, estimated by Equation [2] A clear correlation is seen between the actual and projected densities 
over several orders of magnitude. 



Figure [4] shows a version of the (|Abraham et al.lll996h 
morphology classification system compared to visual 
morphologies. Visual morphologies were provided by one 
of us [RSE] for a complete sample of ~2000 F814W<22.5 
galaxies in the inner, dual band coverage region using the 
precepts discussed in lBundv et a l. (2005). The scatter is 
reasonably large due to the effects of band shifting and 
surface brightens dimming on the asymmetry coefficient 
and visual classification. The difference between Spiral 
and Irregular galaxies is not important for the present in- 
vestigation of the T-S relation; we therefore use only the 
Gini parameter to classify early-type galaxies. We chose 
a cut at logio(Gini)< —0.35 as the dividing line between 
early and late type galaxies because it ga ve a similar frac- 
tion o f early and late type galaxies as the I Abraham et aTl 
(1996) system at 0.2 < z < 0.4. Figure [5] shows a region 
of Figure [4] near the division line between early and late 
type galaxies. The axes are identical, but the symbols 
have been replaced with color images of the actual galax- 
ies. Our specified cut in Gini coefficient clearly separates 
early and late type galaxies. 

4. DENSITY ESTIMATE 

We adopted the N* h neighbor projected density esti- 
mate introduc ed by a number of previous investigators 
(lDresslerlll980HDressler et al.lll997tlPostman et al.ll2005t 
SmitheLaL|l2005). The projected density in the vicin- 
ity of each galaxy is estimated from the distance to the 
N th neighbor. This distance defines the radius of a circle 
whose area is used to estimate the surface density as: 



N 



(1) 



To be consistent with previous studies, we chose to use 
the 10 th nearest neighbor. A median background density 
is then subtracted to correct for line of sight superposi- 
tion due to uncertainties in the distances (ie. redshift). 
This estimator is optimal when line of sight distances 
are uncertain since only projected distances need to be 
estimated. 

Photometric redshifts improve the nearest neighbor 
method in several key ways. The rest-frame absolute 
luminosity of galaxies can be accurately estimated so 



similar galaxies (from the same part of the luminosity 
function) can be selected at all redshifts. The nearest 
neighbor counting can be done in redshift bins rather 
than along the entire line of sight, reducing the back- 
ground contamination and allowing multiple structures 
to be discriminated in the same field. In addition, the 
reduction in background galaxy counts enables one to 
probe to mu ch lower pro j ected densities. 

Following ISmith et al.l (|2005f ) we adopt a luminosity 
cut of My < —21.2 at z = 1 and allow for one magni- 
tude of passive evolution between z = 1 and the present. 
This should select a similar mix of galaxies at all red- 
shifts. Furthermore, maintaining a consistent magnitude 
cut between studies is very important because the galaxy 
densi ty and morphologic al mix may change with magni- 
tude ((Benson et al.ll2001h . 

With this magnitude cut and fading due to band shift- 
ing our F814W data are sufficiently deep at z < 1.4; 
however, a redshift slice of ±3cr z is required around each 
galaxy to ensure all objects in a structure are measured. 
As a result, we are only able to accurately measure den- 
sities for z < 1.3. 

4.1. Accuracy of the Density Estimator 

The accuracy of the photometric redshifts determines 
the minimum density to which the N th nearest neighbor 
method will work. This minimum density can be esti- 
mated by considering a structure of density S s embed- 
ded in a random background of density PBkg, where pskg 
is the number density of galaxies per Mpc 2 per redshift 
interval. The projected density of background sources 
for a slice of thickness Az along the line of sight is then 
Efltj = pBkgAz. Assuming that Az is large enough to 
include all galaxies within the structure, the fraction of 
galaxies actually in the structure is given by: 



Real 



S s + T,Bkg pBkgAz + S s 



(2) 



Equation [2] can also be inverted to define a minimum 
density above which a certain fraction of the galaxies 
will be members of a given structure. We do not need to 
consider Poisson error because £ is determined with the 
same number of galaxies at all densities. 
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In the present data, the largest redshift error is at 
z = 1.3, where a z = 0.065 and E Bfcg = 3Mpc 2 . This 
means 50% of galaxies will be assigned to the correct 
structures at E s = 3Mpc 2 and 77% will be correctly as- 
signed at E s = lOMpc 2 . At z = 0.3, er z = 0.036 and 
^Bkg = 4Mpc 2 , so 42% of galaxies will be assigned to 
the correct structures at E s = 3Mpc 2 and 71% will be 
correctly assigned at E s = lOMpc 2 . 

4.2. Relation of E with Volume Density 

Ideally, one would measure volume densities rather 
than projected densities. Unfortunately, the line of sight 
error from photomet ric redshifts make v olume densities 
difficult to measure (jCooper et al.ll2005h . To first order, 
the relation between the projected mean density, E, and 
volume density is: 

p = B = x [Kl? /2 = 0.42E 1 ' 5 (3) 
H 4r V 167V V ' 

if we assume spherical symmetry. 

To see how well this relation holds we created a sim- 
ulated galaxy catalog with structures of known density. 
The structures had a gaussian density profile and peak 
densities ranging from 1 to 300 galaxies per Mpc 3 , sim- 
ilar to the observed range in COSMOS. The results are 
shown in Figure UK- There is a strong correlation be- 
tween the input and recovered densities for all but the 
sparsest regions (below E ~ 3 galaxies per Mpc 2 ), where 
our assumption of spherical symmetry breaks down. Fur- 
thermore, the same assumption appears to overestimate 
the true density by a factor of ~ 2 — 3. This is also due to 
our assumption of spherical symmetry, which will under- 
estimate the true volume. This offset is only applicable 
if projected densities used in this paper are converted 
to real space densities. No offset is observed between 
projected densities measured with and without redshift 
error. 

In addition to these effects, objects in low density re- 
gions projected in front of or behind a dense region are 
scattered to higher densities; however, the fraction of 
these objects is less than 2% of the total, and can be 
neglected. Figure [6Jd shows the effect of redshift error on 
the density analysis. The scatter is larger than Figure 
due to the line of sight errors, however the real space 
density is clearly r ecoverable. 

Similar tests by ICooper et alj (|2005f) on mock galaxy 
catalogs from CDM simula tions agree with ou r results 
(see Table 2 and Figure 1 in lCooper et all (|2005D ). How- 
ever, the average density in the COSMOS data is signifi- 
cantly higher than those in the Cooper et al. simulations, 
which leads Cooper et al. to conclude that the accuracy 
of photometric redshifts is not sufficient for the majority 
of galaxies. A density of E = 10 galaxies per Mpc 2 cor- 
responds to logio(Ds) — —0.16 in Cooper et al., above 
which there are very few galaxies in their simulation, but 
a significant number of galaxies in the COSMOS data. 

To understand these differences we applied our estima- 
tor to mock galaxy catalogs prov ided by the Millennium 
simulation (Springel et al. 2005). These are the same 
simulation used by Cooper et al., but with dimensions 
and limiting magnitudes for the COSMOS survey. As 
found by Cooper et al. densities above E = 3 galax- 
ies per Mpc 2 are successfully measured with photometric 
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Fig. 7. — The number of galaxies per unit area at each density 
are shown for COSMOS and mock catal ogs from two versions of the 
Millennium Simulation (Springel et al. 2005). The vertical dashed 
line indicates the minimum measurable density, S = 3, estimated 
by Equation [2] The v2.0 mock catalog populates the dark matter 
haloes in a similar way to those used in lCooper et all (2005) while 
the v3.0 mock catalog follows galaxy orbits in detail. The selection 
function, redshift error, and area of the simulations are identical to 
the COSMOS data. Notice the distribution of densities in the v2.0 
mock catalog is significantly lower than the actual data. Error bars 
indicate the measurement Poisson error except for the v3.0 mock 
catalog which also includse the expected range of cosmic variance. 

redshifts. However, the number of galaxies at densities 
greater than 3 galaxies per Mpc 2 in the Millennium sim- 
ulation is considerably lower than those in the COSMOS 
data. Figure [7] shows the distribution of densities re- 
covered from the Millennium simulation and from the 
COSMOS survey. At densities above E > 3 there are 
a factor of ~ 10 more galaxies in the COSMOS survey 
than the v2.0 simulations. 

The discrepancy in density distribution is due to the 
way galaxies are distributed within dark matter haloes. 
In earlier versions of the Millennium simulation galaxy 
orbits were not followed in detail after dark matter haloes 
merged (M. Kitzbichler & S. White private communica- 
tion). This resulted in systematically fewer galaxies in 
high density regions. This has been corrected in the lat- 
est versions (v3.0 and newer) of the simulation, improv- 
ing the agreement, but still underestimating the number 
of gal axies in high dens i ty reg ions (See Figure [7]). How- 
ever, iMcCracken et al.l (|2007l ) find a higher amplitude 
in both the overall correlation function and the correla- 
tion function on small scales than predicted by the mock 
catalogs, which indicates the mock catalogs still tend to 
under-populate dense regions. 

5. COMPARISON WITH LITERATURE DATA 

It is important to maintain a consistent magnitude 
limit when studying the morphology density relation be- 
cause the galaxy densi ty and morphologi cal mix may 
ch ange with magn i tude dBenson et al.ll2001l) (see Section 
4). iDresslerl (119801) ^[Dressier et al I (|1997D ; lPostman et al.l 
(2005) and iTreu et alj (|2003l ) use a magnitude limi t 
more than a magn itude fainter than iGoto et al.l (|2003f) : 
ISmith et all (|2005f ) and our work (see Figure HJ). To 
quantify the effects of the different limiting magnitudes 
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TABLE 2 

Parameters of Previous T-E Studies 



Reference 



Redshift 


H 




Qy 


Magnitude Limit 


Correction to 


Correction 


Range 


(km s _1 Mpc -1 ) 






M v < 1 


Density 2 


F £+S0 


0.011-0.066 


50 


1.0 


0.0 


-19.75 


-0.33 


0.0 


0.05-0.1 


75 


0.3 


0.7 


-20.3 3 


0.0 


0.0 


0.37-0.56 


50 


1.0 


0.0 


-20.0 


-0.47 


-0.045 


0.4 


65 


0.3 


0.7 


-18.7 


-0.38 


-0.056 


0.78-1.27 


65 


0.3 


0.7 


-21. 2 4 


-0.12 


0.0 


0.4-1.27 


70 


0.3 


0.7 


-20.07 4 


-0.32 


-0.045 



Dressier (1980) 
Goto et al. (20031 
Dressier et al. (19971 
Treu et al"72003T " 
Smith et al. (20051 
Postman et al. (20051 



1 In a Ho — 75, £T m — 0.3, Hy — 0.7 cosmology 2 In units of logio(£). The correction includes conversion to a Ho — 75, 
S^m — 0.3, f2v — 0.7 cosmology and an offset for the different magnitude limits. 3 The actual limit is M r < —20.5, 
which corresponds to My < —20.3 at the median galaxy color. 4 The values given are for z — 1, ISmith et al.l J2005H 
allows for one magnitude of passive evolution between 2—1 and z = 0, while IPostman et all 1 120051) allows for 0.8 
magnitude of passive evolution. 




I 

0.2 



I 

0.4 



COSMOS 

Dressier et al. 1980, 1997 
Goto et al. 2003 
Postman et al. 2005 
Smith et al. 2004 
Treu et al. 



1 1.2 



The density measured with both fainter magnitude 
limits is 0.26 dex higher than that measured with My < 
—21.2 at densities above S > 3 galaxies per Mpc 2 . This 
offset appears to vary with density at S < 3 galaxies per 
Mpc 2 . However, these densities are not reliable (see Sec- 
tion 4.1) and the trend is not apparent at higher density. 
No trend with redshift was observed in either magnitude 
bin or at any density. 

The early-type fraction is 0.045 lower using a limit of 
My < -20.05, and 0.056 lower with a limit of My < 
— 19.05 than with a limit of My < —21.2. No significant 
trends in these offset are observed with density between 
3 < £ < 100 galaxies per Mpc 2 or redshift between 0.2 < 
z < 0.8. However, we do not probe the highest densities 
and redshifts used in the literature. 

A summary of the literature data and the correction 
factors applied is given in Table [2] 



Fig. 8. — The limiting magnitude for various studies of the T-E 
relation is shown along with the limit used in this sample (Solid 
Line). The Goto et al. (2003) point has been converted from M T 
to My assuming the median rest frame color obse rved in the COS - 
MOS data between 0.2 < z < 0.4 ((V - r) = 2). IDresslerl lli980l 1; 
Dressier et al. (1997) and Postman et al. (2005) use a limiting mag- 
nitude 1 magnitude fainter than our study. The dashed line indi- 
cates a limiting magnitude of My < —19.05 with one magnitude 
of passive evolution. 



in the literature we began by constructing morpholog- 
ically selected luminosity functions uncorrected for in- 
completeness at the faint end. These are shown for five 
redshift bins in Figure The solid lines mark our mag- 
nitude limit, while the dashed line indicates the fa inter 
limit used by [D ressier (1980); D ressier et al.l (|1997l ) and 
IPostman etaLI (|2005l ). A dotted line in the 2 < z < 0.4 
bin indicates the limit used by iTreu et all ([2003) . The 
overall fraction of late type galaxies is clearly higher at 
fainter magnitudes. In addition our data becomes incom- 
plete at z > 0.8 for the fainter magnitude limits. 

We estimate corrections for the differing magnitude 
limits by analyzing our data with three magnitude lim- 



= 1: My 



< 



-21.2, My < -20.05, and 



its at z 

My < —19.05. All three analysis include 1 magnitude of 
passive evolution between z = 1 and the present. The 
My < -20.05 limit corres ponds to lDressler et al.l (|1997l ) 
an d IPostman et alj (120051 ) . while My < -19.05 matches 
the lTreu et al l (|2003l ) work at z = 0.4. 



6. RESULTS 

Our measurements of the T-S r elation 
with those taken from the literature (iDre sslc 



Dressier et al.l 



Postman et al 



along 
_3[l980; 

19971: iGoto et~all[200l ISmith et all 12005 : 
2005), are shown in Figure [10] for five red- 



shift bins. The densities for the literature points have 
been converted to a fi„ = 0.7, O m = 0.3, H a = 75 cos- 
mology and corrected for differences in the limiting mag- 
nitude. After applying these corrections, our data is con- 
sistent with the earlier studies. The COSMOS data sig- 
nificantly improves the precision of the early-type mor- 
phological fraction, and increases the redshift and den- 
sity resolution, compared to earlier studies. We find that 
the T-E relation was already in place at z > 1, but dif- 
fers from the local re l ation. As seen in previous studies 
([Postman et al.ll2005l : ISmith et al1l2005h . the early- type 
fraction is smaller and increases more gradually with den- 
sity at z = 1 than the at z = 0. 

The evolution of the T-E relation with redshift is en- 
capsulated in Figure 1111 This figure shows the fraction 
of early-type galaxies as a function of time for four differ- 
ent density bins — growth in the early-type fraction with 
look-back time is seen at all densities with more rapid 
growth at higher densities. However, the majority of the 
evolution occurs at z < 0.4 for densities below ~ 100 
galaxies per Mpc 2 . Assuming the growth rate of the 
early-type fraction is smooth, a line can be fit to the 
data at each density bin in Figure [TT] The slope of this 
line yields the growth rate of the early-type fraction at 



0.1 - 



0.01 =■ 



0.001 e- 



0.0001 - 



0.2<z<0.4 



0.1 = 



0.01 - 



0.001 =■ 



0.0001 - 



0.1 



a o.oi 
■a 

9) 



b 
p. 

* 



0.001 = 



0.1 



0.01 



0.001 



0.0001 =■ 



0.1 - 



0.01 ; 



0.001 ; 



0.0001 : 



-25 -24 -23 -22 -21 -20 -19 

NL, 



-18 











E 0.4<z<0.6 


* 5 : 





































-25 -24 -23 -22 -21 -20 -19 

NL, 



-18 



: 0.6<z<0.8 


i 


=©= =0= 




* ^ * 1 


— a 


m 



















-25 -24 -23 -22 -21 -20 -19 

NL, 



-18 









= 0.8<z<1.0 






_ -6 


-e- -©- 


-e- 


* * * 




=©= — 








_ 36: 






BC 












- ^ 













-25 -24 -23 -22 -21 -20 -19 -18 

NL, 



1.0<z<1.2 



Early Type 
Late Type 



1 

0.8 
0.6 
0.4 
0.2 
0.8 
0.6 
0.4 
0.2 
0.8 
i 0.6 h 



0.4 
0.2 
0.8 
0.6 
0.4 
0.2 
0.8 
0.6 
0.4 
0.2 



Elliptical Fraction 

By density and redshift 




0.6 z<0.8 




1 2 



O Cosmos 

A Dressier et. al. 1980, 1997 
x Goto et. al. 2003 
O Postman et. al. 2005 
O Smith et. al. 2004 
Treu et al. 2004 



-25 -24 -23 -22 -21 -20 -19 -18 

NL, 



Fig. 10. — The morphology density relation is show for five 
redshift intervals, 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8, 
0.8 < z < 1.0, and 1.0 < z < 1.2, along with results from similar 
studies. The literature values have been converted to a cosmol- 
ogy with Q v = 0.7, Q m = 0.3, and H = 75 and corrected for 
differences in the limiting magnitude of the samples. Note the in- 
crease in Elliptical fraction with density and dec re asing redshift. 
These results agree with those of IDressl erl dl980 | ) : [Dressier et al.l 
1119971 ); [Goto et al.l IL2003T) ; IPostman et al.l 112003 ). The lowest den- 
sity point from ISmith et alj (2005J) is discarded due to their large 
redshift error (see Section 4.1). 



Fig. 9. — The luminosity function of galaxies in COSMOS split 
by morphological type is shown for five redshift bins. These are 
uncorrected for incompleteness. The solid lines mark our absolute 
magnitude limi t, the dashed lin e indic ates t he fainter limit used b y 
IDresslerl fl980T) ; IDressler etaLl ( fl997l> and IPostman et~aH (|2005h . 
and the dotted l ine in the top panel indicates the limit used by 
ITreu et al.l II2003Q . Notice the larger fraction of late type galaxies at 
fainter magnitudes and the fact that our data becomes incomplete 
at z > 0.8 for the fainter absolute magnitude limit. 
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Lookback Time (Gyr) 



Fig. 11. — Evolution of the elliptical fraction with redshift 
is shown for four densities. Notice the increased rate of ellipti- 
cal formation with increased density and the lack of evolution at 
z > 0.4 i n the lowest two density bins. Open p o ints are a compila - 
tion fromlDresslerl fl98Ch;IDressler et al.l lfl997h ; ITreu et alj (2003); 
IGoto et all H2003I 1: [Postman et al.l ^200511 and ISmith etTaO (|2005h , 
solid points are those from this study, the lines are best fits to the 
data. 



0] 











O All Data 
X z>0.4 data 

Best Fit Slope = 0.29+0.02 

No Evolulion 










T C 


) 













\.\ 1 10 100 1000 10000 

Density (Galaxies Mpc ) 



Fig. 12. — The growth rate of the E+S0 fraction as a function of 
density is plotted. The thick solid line is the best fit power law to 
the data with a slope of 0.29±0.02, while the dashed line shows the 
expected trend for no change in evolution rate with density. The 
red-points indicate the measured evolution if only data at z > 0.4 is 
used. No evolution is observed at densities below E < 100 galaxies 
per Mpc 2 at these high redshifts. These rates are derived from all 
data p rese nted in Figure 1101 not the selected densities plotted in 

Figure rrn 

any given density, which are plotted in Figure 1121 

Figure [12] shows a clear trend for more rapid early-type 
(E+S0) production at higher densities. The data are well 
fit by a power law with a slope of 0.29 ± 0.02. However, 
without the z < 0.4 points the growth rate drops off 
quickly at densities below ~ 100 galaxies per Mpc 2 . It 
should be noted that points above £ > 100 galaxies per 
Mpc 2 come entirely from cluster surveys in the literature. 
However, it is unlikely that the observed evol ution is due 
to sy s tematic effects because all three studies dTreu et al.l 
l2003t iSmith et al.l [20051 : iPostman et all ;2005) at these 
redshifts and densities use the same visual morphological 



O Morphological Type 



□ Spectral Type 
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Fig. 13. — The growth rate in the non-star-forming (passive) 
galaxy fraction (blue) is shown along with the growth rate in the 
early-type morphological fraction (black) as a function of density. 
The measured rates assume the growth rate is linear with time. 
The thick solid line is the linear growth rate expected by the best 
fit closed box model (Equations [5] and [6j at the weighted mean 
lookback time of the COSMOS data (6.5Gyr). The best fit model 
has Y/ = 1.07S 08 and t = -0.99 * log 10 (T,) + 1.77Gyr. Notice the 
growth in the non-star-forming galaxy fraction is well fit by the 
closed box model, however no reasonable set of parameters can fit 
the growth in the morphological fraction unless more early-types 
are intrinsically formed in dense regions. 



classification system, and the tre nd to stronger evo lution 
at higher densi t y is se en by both ISmith et al.l (|2005f ) and 
IPostman et al.l (|2005l ) independently. 

The growth in the fraction of non-star forming (pas- 
sive) galaxies is shown in Figure [13] along with the data 
from Figure 1121 The spectral types are drawn from a 
similar analysis in IScoville et al.l (|2007al) . however, the 
density measure and magnitude l imit are the same one s 
used in this study, not those from IScoville et all (|2007al ). 
Note the difference in the evolution measured from star 
formation and morphology. 

Many bulge dominated but elongated objects with 
higher asymmetry (seen near the top of Figure 2]) are 
selected by both the Abraham et al. and our early-type 
selection. These elongated objects are likely edge on SO 
galaxies. The Gemini-Deep-Deep-Survey (GDDS) finds 
that many of these elongated ob jects have residual star 
formation (| Abraham et alj |2004). Although such object 
with active star formation activity do not strictly meet 
the criteria of early-type galaxies, they are bulge domi- 
nated, and hence more dynamically relaxed than galax- 
ies with smaller Gini co-efficent. These objects are likely 
progenitors of SO's or galaxies with early-type morphol- 
ogy undergoing a burst of star formation. Furthermore, 
since our Gini parameter is invariant with redshift, the 
same fraction of the galaxy population will be selected 
at all redshifts, implying that any contamination will not 
affect the observed evolution. 

However, it is important to remember there are multi- 
ple sub-classes of early-type g alaxies, and each sub-c lass 
could be evolving differently. IPostman et al.l (|2005f ) in- 
vestigate this effect by analyzing the ellipticity distribu- 
tion of cluster early-type galaxies and find no evolution 
in the ellipticity distribution between z ~ 1 and z ~ 0.6. 
A similar analysis of our data for S > 10 galaxies per 



10 



Mpc 2 also finds no evolution in the ellipticity distribu- 
tion of early-types. A Kolmogorov-Smirnov test finds 
the distribution of cllipticities at 0.2 < z < 0.4 is consis- 
tent with those at 0.8 < z < 1.2 with 96.5% certainty. 
Therefore, any differential evolution in the E and SO pop- 
ulation is not seen in the overall ellipticity distribution 
of early- type galaxies. 

7. DISCUSSION 

Whether the morphology-density (T-E) relation is an 
intrinsic property of cluster galaxies or a result of envi- 
ronmental influence has been a matter of debate for some 
time. In contemporary simulations of galaxy formation 
the majority of galaxies begin as dynamically cold, gas 
rich, star forming disks. They are subsequently trans- 
formed into early-type systems through interactions with 
other galaxies or by losing or exhausting the gas required 
to form stars. In this context there are two hypothesis ex- 
plaining the T-E relation: dense regions formed earlier 
than sparse ones, and dens regions have stronger and 
more frequent interactions. A third hypothesis is that 
dense regions intrinsically form a higher fraction early- 
types, and that no transformation is needed. 

The relative importance of these hypothesis can be sep- 
arated with a simple closed box model. If galaxy clusters 
were closed systems one would expect the growth in the 
early-type fraction to slow with time as the fraction of 
early-type galaxies increases. This occurs because there 
are fewer late type galaxies left to transform as the early- 
type fraction approaches 100%. Formally, the growth 
rate is given by: 



dNE+S0 =(1-^±^)JVt,(E) 



dt 



N 



(4) 



where the rate of change in the number density of early- 
type galaxies Ne+so is given by the total number density 
of galaxies N, and the conversion rate of late to early- 
type galaxies per unit time as a function of density, 77(E). 
Integrating Equation [3] with respect to time we get: 



Fe+so = 1-(1- F„(Ti))e 



-r,(S)(t-t ) 



(5) 



where Fe+so is the early-type fraction at time t, F (E) is 
the primordial fraction of early types at a given density 
and t is the cluster formation time. This equation can 
be differentiated to yield: 



dF, 



E+SO 



dt 



= (1 - F o (E))77(E)e-" (s)(t -* 0) 



(6) 



If we assume 77(E) =constant and F G (E) = the 
observed increase in the production rate of early-type 
galaxies with density suggests cluster galaxies formed 
significantly later than the field galaxies. Since de- 
tailed studies of cluster ellipticals indic ate they are 
older than their coun t erparts in the field (iBower et al 



1992t lEllis et al.lll997t IStanford et alJll998t jLucey et all 
199l[ Ivan D okkum fc Framd 119961: iBender et all 119981 : 



Pahre et all H998: Kels on et al J 120001 ). we can rule out 
differences in formation time as the sole source of the 
T-E relation. A model were 77(E) increases with den- 
sity, but F (Y.) = const, is also ruled out by the data if 
we require that formation time is constant or increasing 
with density. A combination of earlier formation times in 



dense regions, increasing transformation rates with den- 
sity, and an intrinsically higher early type fractions in 
dense regions is required to explain the observed data 
with a closed box model. However, simply allowing for 
galaxy in- fall also explains the observed e volution . 

Dynamical frict ion or "harassment" fSp itzerl 119581 : 
IMoore et all I1998D is the most likely source of an in- 
creased late to early type transformation rate in dense 
regions. This process operates by increasing the orbital 
energy of stars inside individual galaxies through tidal 
interactions. The rate of momentum exchange due to 
dynamical friction increases with density, which natu- 
rally leads to the qbsery e d inc rease in transformation 
rate. Furthermore, iGotol (|2005l ) find differences in the 
velocity dispersion of cluster early and late type galaxies 
which can be explained if the early type galaxies under- 
went more dynamical friction than the late type galax- 
ies, but can not be explained by gas stripping. Finally, 
the amount of energy transfered to a galaxy through 
dynamical friction is proportiona l to the gala xy mass 
(|Spitzerl H951I IMoore et al.l 119981 IGotol l200l . so the 
color-magnitude and fundamental plain relations may 
also be a consequence. 

Intrinsic differences in the galaxy formation process be- 
tween clusters and the field also appear to be important 
in determining galaxy morphology. Using only our data, 
no evolution is measured in the T-E relation, but COS- 
MOS only covers redshifts greater than 0.4 and densities 
below E < 100 galaxies per Mpc 2 . At hig her densities, 
evolut ion is independently observ ed by both lSmith et all 
(2005) and (|Postman et al.ll2005f ). Nevertheless, the lack 
of evolution at low densities suggests the galaxy forma- 
tion process may be intrinsically different in dense re- 
gi ons than sparse. 

iBundv et al.l (|2006l ) also find intrinsic differences in the 
galaxy formation process may be more important than 
environment. They show the star formation properties 
of galaxies are more correlated with galaxy mass than 
environment. Environment is only important when the 
local galaxy density was significantly greater than the 
field density. If a similar statement can be made for 
morphologies, it would explain why evolution is only ob- 
served in the highest density regions. 

However, morphology and star formation appear to be 
affected by different processes. The fraction of galaxies 
with early-type morphologies and the fraction with low 
star-formation rates evolve in different ways with den- 
sity (see Figure |T3j) . These different rates of evolution 
can be explained if the growth in the early type fraction 
were driven by interactions, while the reduction in star 
formation was caused by gas stripping . 

Assuming gas removal due to cluster interactions is 
responsible for truncating the star formation, we can re- 
cycle our closed box model to describe the process. For 
the SFR-E relation F (Y.) — because stars must form 
to make galaxies. If the gas density traces the galaxy 
density, 77(E) becomes the rate at which star-formation 
is truncated and t becomes the formation time for the 
cluster. This simple model accurately reproduces the ob- 
served shape and evolution of the SFR-E relation while 
also predicting reasonable galaxy formation times. How- 
ever, points at much lower density than what we can 
probe with photometric redshifts are required to confirm 
the gas striping hypothesis. 
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Nevertheless, several studies of the local universe also 
suggest the SFR-S relation is driven by gas stripping, 
while the T-S relat ion is due to galaxy interactions. 
iBlanton et al.l (|2005[ ) find color is a better predictor of 
density than morphology, indicating star formation is 
strongly affected by environment while early-type mor- 
phologies are not necessarily the result of environment . 
Further evidence is provided by iQuintero et all (2006), 
who find both morphology and color are correlated with 
distance from cluster centers, but the correlation between 
morphology and star-formation is asymmetric. Specifi- 
cally, for a given star formation rate, the fraction of mor- 
phological early-type galaxies does not change with dis- 
tance from the cluster, but at a given morphological type, 
the average star formation rate in creases with d i stance 
fro m the cluste r . In a ddition, both Dres sier et alj (|1997f ) 
and lTreu et alj (|2003l ) show morphology is more strongly 
correlated with local density than distance to the cluster 
center, indicating the number of neighbors and hence the 
number of interactions is the more important quantity in 
determining morphology. 

Galaxy in-fall may also play a role in shaping the T-S 
relation. It is predicted by CDM models, but it is difficult 
to explain the small scatter in the colors and fundamen- 
tal plane of cluster ellipticals if most early-type galaxies 
were recently transform ed from field galaxies. Neverthe- 
less, iTreu et all (|2003[ ) point out that the dynamics of 
galaxy clusters would erase the T-S relation if the or- 
bits of individual galaxies were not confined to regions 
of nearly constant density. So, if in-fall is responsible 
galaxies are likely accreted smoothly via dynamical fric- 
tion from adjacent regions of already elevated density. 

There is no quantitative study of how harassment or 
gas sta rvation should affect the evolution of the T-S re- 
lation. iBenson et al.1 (|2001f ) attempts to model evolution 
of the T-S relation at low densities in cold dark matter 
(CDM) simulations via halo mergers. The amplitude of 
the predicted evolution is correct at S = 10 galaxies per 
Mpc 2 , and no differential evolution is seen at lower den- 
sities. Unfortunately, at higher densities, where we have 
the bulk of our data, the Benson et al. model is not valid 
because individual galaxy orbits are not tracked once the 
dark matter halos merge, so no comparison can be made. 

Clearly more detailed models are needed, but care 
must be taken when comparing the results. A 
magnitude-morphology relation is predicted by most 
CDM models, so changing the magnitude limit used to 
measure density will change the observed T-S relation 
and its observed evolution. As a result, it is very impor- 
tant to match the absolute magnitude cuts when com- 
paring models and data. Furthermore, the present mea- 
surements do not distinguish between moderate density 
regions on the outskirts of massi ve clusters and the cen- 
ters of moderate mass groups. ITreu et al.l rf2003h find 
this distinction is not important, nevertheless, it remains 
a potential source of systematic uncertainty in our mea- 
surements. 

In considering Figure lll| several limitations of our 
present data become evident. Our measurement of the 
early-type growth rate strong ly dep end on the local mea- 
surements from iGoto et al.l (|200a ) and iDresslerl (|1980l ) 
which use a different morphological classification scheme. 
Furthermore, the large gap in data between 0.1 < z < 
0.3, which is nearly a quarter of the cosmic time we 



probe, limits our ability to measure how the growth in 
the early-type fraction changes with time. COSMOS has 
too small an area and the SDSS has neither the depth 
nor the spatial resolution to probe the T-S relation in 
this gap. To probe a volume at 0.1 < z < 0.3 compara- 
ble to COSMOS at z = 0.5, a survey would have to cover 
« 15 square degrees with 0.6" seeing in one band. 

8. CONCLUSIONS 

We have developed a technique for measuring the 
early- type galaxy fraction with density and redshift. The 
Gini parameter is found to reliably select early-type 
galaxies between 0.3 < z < 1.2 using only the ACS 
F814W filter. This selection is free from systematic ef- 
fects in red-shift (band shifting) or surface brightness. 

We find densities are measurable with photometric red- 
shifts if projected densities are used. Nevertheless, for 
photometric redshift accuracies of A z /(1 + z) = 0.03 
projected densities below S = 3 galaxies per Mpc 2 are 
difficult to measure. The observed number of galaxies at 
a S > 10 galaxies per Mpc 2 is significantly higher than 
that found in earlier CDM galaxy simulations because 
these simulations did not follow the orbits of galaxies 
within dark matter haloes. The latest simulations now 
follow galaxies, improving the agreement, but still un- 
derestimating the number of galaxies in dense regions. 

Using these techniques we measure the evolution of 
the T-S relation and the growth rate of the early-type 
fraction with cosmic time. We find the growth rate of 
the early-type fraction is increasing with density, which 
can not be explained by a closed box model with early 
formation times. We conclude some density dependent 
process combined with galaxy in fall is responsible for the 
observed relation and evolution, with dynamical friction 
and harassment being the most likely mechanisms. 

The SFR-S relation appears to result from different 
processes than the T-S relation. In particular, SFR-S re- 
lation evolves differently from the T-S relation and there 
does not appear to be a direct relationship between the 
two. Gas stripping is the most likely source of the SFR-S 
relation. 

At S = 10 galaxies per Mpc 2 , the rate of dark halo in- 
teractions in CDM models predict the correct amplitude 
of growth in the early-type fraction, but these models are 
not valid at higher densities. No current model attempts 
to predict our observed evolution at S > 10, so we can 
only speculate as to what physical mechanisms must be 
at play. Furthermore, the current generation of CDM 
models appear to under- estimate the number of galaxies 
in high density regions. 
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